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Abstract 



In this paper we study properties of Lax and Transfer matrices associated with quantum inte- 
igrable systems. Our point of view stems from the fact that their elements satisfy special commutation 
^ ^properties, considered by Yu. I. Manin some twenty years ago at the beginning of Quantum Group 

I ^Theory. They are the commutation properties of matrix elements of linear homomorphisms be- 

<^ 'tween polynomial rings; more explicitly they read: 1) elements in the same column commute; 2) 
.commutators of the cross terms are equal: [Mij,Mki] = [Mkj,Mii\ (e.g. [Mii,M22] = [M21, M12]). 
! The main aim of this paper is twofold: on the one hand we observe and prove that such matri- 
^ 'pes (which we call Manin matrices for short) behave almost as well as matrices with commutative 
^ 'elements. Namely theorems of linear algebra (e.g., a natural definition of the determinant, the 
'""''Cayley-Hamilton theorem, the Newton identities and so on and so forth) have a straightforward 
psj .counterpart in the case of Manin matrices. 

^ '. On the other hand, we remark that such matrices are somewhat ubiquitous in the theory of 
^ Iquantum integrability. For instance, Manin matrices (and their q-analogs) include matrices satisfying 
^ the Yang-Baxter relation "RTT=TTR" and the so-called Cartier-Foata matrices. Also, they enter 
• 0?alalaev's' remarkable formulas: det{dz — Lcaudiniz)), det{l — e~^'TYangian{z)) for the "quantum 
^ spectral curve" , and appear in separation of variables problem and Capelli identities. 

! We show that theorems of linear algebra, after being established for such matrices, have various 
lapplication^o quantum integrable systems and Lie algebras, e.g in the construction of new generators 
in Z (Ucritigln)) (and, in general, in the construction of quantum conservation laws), in the Knizhnik- 
^ -Zamolodchikov equation, and in the problem of Wick ordering. 

. We also discuss applications to the separation of variables problem, new Capelli identities and 
'the Langlands correspondence. 
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1 Introduction and Summary of the Results 



It is well-known that matrices with generically noncommutative elements do not admit a natural 
construction of the determinant and basic theorems of the linear algebra fail to hold true. On 
the other hand, matrices with noncommutative entries play a basic role in the theory of quantum 
integrability (see, e.g., [FTS79]), in Manin's theory of "noncommutative symmetries" |Man88j . and 



so on and so forth. It is fair to say that recently D. Talalaev [Ta04] made a kind of breakthrough 
in quantum integrability, defining the "quantum characteristic polynomial" or "quantum spectral 
curve" det{dz — L{z)) for Lax matrices satisfying rational i?-matrix commutation relations (e.g., 
Gaudin systems jl. 

The first and basic observation of the present paper is the following: the quantum Lax matrix 
[dz — Lgi^^caudin{z)) entering Talalaev's formula, as well as a suitable modification of the transfer 
matrix of the (XXX) Heisenberg chain, matches the simplest case of Manin's considerations (i.e. 
are Manin matrices in the sense specified by the Definition [1] below). Further we prove that many 
results of commutative linear algebra can be applied with minor modifications in the case of "Manin 
matrices" and derive applications. 

We will consider the simplest case of those considered by Manin, namely - in the present paper 
- we will restrict ourselves to the case of commutators, and not of (super)-g-commutators, etc. Let 
us mention that Manin matrices are defined, roughly speaking, by half of the relations of the 
corresponding quantum group Funq{GL{n)) and taking q = 1 (see remark [2] page [6]) . 

Definition 1 Let M be a n x n' matrix with elements Mij in a noncommutative ring TZ. We will 
call M a Manin matrix if the following two conditions hold: 

1. Elements in the same column commute between themselves. 

2. Commutators of cross terms 0/2x2 submatrices of M are equal: 

[M,„ Mki] = [Mfc,-, Mai V2, J, fc, / e. g. [Mn, M22] = [M21, M12]. (1) 

A more intrinsic definition of Manin matrices via coactions on polynomial and Grassmann algebras 
will be recalled in Proposition [1] page [7] below. (Roughly speaking variables Xi = Ylj ^ij^j commute 
among themselves if and only if M is a Manin matrix, where Xj are commuting variables, also 
commuting with elements of M). 

As we shall see in Section |3], the following properties hold: 

• dz — Lgi^_Gaudin{z) is a Mauiu matrix. Here Lgi^_Gaudin{z) is the Lax matrix for the Lie algebra 
gln[t], i-e. a matrix valued-generating function for generators of gln[t], (in physics Lgi^_Gaudin{z) 
is the Lax matrix for the Gaudin integrable system). (Section 13.11) . 

• e~^''Tgi^_Yangian{z) IS a Mauiu matrix. Here Tgi^_Yangian{z) is the Lax (or "transfer") matrix 
for the Yangian algebra Y{gln), i.e. a matrix valued-generating function for the generators 
of Y{gln), (in physics Tgi^_Yangian{z) is the Lax matrix for the integrable systems with the 
Yangian symmetry: the Heisenberg XXX spin chain, the Toda system, etc.) (Section 13.21) . 

Furthermore, we shall see that Manin matrices enter other topics in quantum integrability, e.g, the 
quantum separation of variables theory (section [6]) and Capelli identities (section 14.3. ip . 
Among the basic statements of linear algebra we establish (in a form identical to the commutative 
case) for Manin matrices, let us mention: 



See [CT06][ [MTV06][ [MTV05] for some applications to the Bcthc ansatz, separation of variables, the Langlands 
correspondence, the Capelli identities, real algebraic geometry and related topics. 

^Remark that Lgi^-Gaudin{z) is not a Manin matrix; to insert dz is an insight due to D. Talalaev; the operator 
{dz — Lgi^-Gaudin{z)) is also related to the Knizhnik-Zamolodchikov equation see section [4. 1.11 



3 



• Section 14.21 The inverse of a Manin matrix M is again Manin; furthermore det{M ^) = 
{det M)-i; 

• Section I4.3t the formula for the determinant of block matrices: 

det(^^ D ) " rfet(A)det(D - CA~^B) = det{D)det{A - BD-^C); 

Also, we show that D — CA~^B, A — BD^^C are Manin matrices; 

• Section [5?T} the Cayley-Hamilton theorem: det{t — M)\t=M = 0; 

• section [521 the Newton identities between TrM^ and coefficients of det{t + M)). 

This extends some results previously obtained in the literature: Manin has defined the determinant, 
proved Cramer's inversion rule, Laplace formulas, as well as Plucker identities. In [ GLZ03 the 
MacMahon-Wronski formula was proved; [Ko07BJKo07A contains Sylvester's identity and the 
Jacobi ratio's theorem, along with partial results on an inverse matrix and block matricesj^l One of 
our point is to present applications to the theory of Quantum integrable systems. 



Applications to integrability and Lie theory. We derive various applications of the above- 
mentioned linear algebra statements to quantum integrable systems, to the Knizhnik-Zamolodchikov 
equation, as well as Yangians and Lie algebras. Namely: 



Section 14.1.11 a new proof of the correspondence ( CT04 ) between solutions of the Knizhnik- 
Zamolodchikov equation and solutions Q{z) of det{dz — L(z))Q{z) = 0. The last equation plays 
a key role in the solution of the Gaudin type integrable systems. In the s/2 case it is known in 
the literature as "Baxter's T-Q equation". 

Section 14.2.11 the observation that the inverse of a Manin matrix is again a Manin matrix 
yields, in a very simple way, some recent result on a general construction about quantum 
integrable systems ( [E ROl] , [BT02 ] ), that, in a nutshell, says that (quantum)separability implies 
commutativity of the quantum Hamiltonians. 



Section 14.3.11 a new proof of the generalized Capelli identities MTV06 for the Lie algebra 
gl [t] (or " the Gaudin integrable system" ) . 

Sections 15.2. 1^ 15.2.21 a construction of new explicit generators (: TrL[^l(z) :) in the center of 
Uc=critigin[t,'t~^], c) and {TrL^''^{z)) in the commutative Bethe subalgebra in U{gln[t\); in the 
physics language {TrL^''^{z)) are quantum conservation laws for the quantum Gaudin system. 
Here L^''^{z) are "quantum powers" defined in [CT06b] , that "cures" some bad properties (i.e.. 



failure of producing quantum commuting operators) of the standard powers L^{z). Similar 
results for the Yangian case, and also Newton and the Cayley-Hamilton identities for these 
generators. 



^These authors actually considered more general classes of matrices 

^The basic definitions on "normal ordering" : ... :, "critical level", Uc=crit{gln) = Uc=crit{gln[t,t~^], c), and 
("Bethe") subalgebra in U{gln[t])) arc recalled in the Appendix (section [HH]) 

''E.g., [TrL'^gi^_Q^^^^^{z),TrL^gi^ -Oaudzni^)] ¥= ^[CRyT04]j this is expected to be a general phenomenon in the 
quantum case, opposite to the classical one. 
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• Section 15.3.11 the construction of a further set of exphcit generators (: TrS^L{z) :) and 
{TrS^ L[z)) in the center of ?7c=crjt(fl'^n[^, ^~^], c), and, respectively, in the commutative Bethe 
subalgebra in U{gln\t\)] in physical language, (TrS^L{z)) are quantum conservation laws for 
the quantum Gaudin, Heisenberg, etc. integrable systems. We notice that some results on 
the Cayley-Hamilton, Newton, MacMahon-Wronski relations in U{gln) by T. Umeda, M. Itoh 
p98]jU03], .[Ito98], are implied by our constructions. These issues will be discussed elsewhere. 

Finally, we address the problem of Sklyanin's separation of variables problem; we present a conjectural 
construction of the quantum separated variables for Yangian type systems, as well as some remarks 
concerning the Gaudin case. 

In the Appendix, we first collect definitions about the center of Uc=crit{g^n[t,t~^],c), that is, the 
commutative Bethe subalgebra in U{gl[t]) (Section 18.11) . and on the relations of the present work 
with the geometric Langlands program (Section 18.21) . In Sections 18.31 and [6l we give a new look and 
generalizations of the Capelli type identities [Ok96] , [MTV06] , which appears to be related with the 
quantization of the Gaudin and Heisenberg's XXX systems via the Wick ordering prescription. 

Our starting point is the following remarkable construction by D.Talalaev of the "quantum 
spectral curve" or the "quantum characteristic polynomial". It solves the longstanding problems 
of the explicit efficient construction of the center of Uc=crit{g^n[t,t~^],c), the commutative Bethe 
subalgebra in U{gl[t]) |^ and explicitly produces a complete set of quantum integrals of motion 
for the Gaudin system. It also has many other applications (to the Bethe ansatz, separation of 
variables, the Langlands correspondence, the Capelli identities, real algebraic geometry etc). The 
right setup where the Talalaev's formula fits are the ideas of E. K. Sklyanin on quantum separation 



of variables (see, e.g., the surveys ![Sk92], [Sk95] ) 



Theorem [Ta04] Let L{z) be the Lax matrix of the gln-Gaudin model, and consider the following 
differential operator in the variable z (the Talalaev's quantum characteristic polynomial or quantum 
spectral curve): 

def'''^^''{d,-Lgi^^Gaudin{z))= QHi{z)di, (2) 

i=0,...,n 

Then: 

Vi,j G 0, ...,n, andu eC,v eC [QHi{z)l=u,QHj{z)l=.y] = 0. 

So taking QHi{z) for different i, z = u ^ C (or their residues at poles, or other "spectral invariant" 
of the Lax matrix dz — Lgi^^Gaudin{z)) one obtains a full set of generators of quantum mutually 
commuting conserved integrals of motion. 

Actually, the theorem holds for all Lax matrices of "(yf/^-Gaudin type (see definition 0, section [XT!) . 
For the Yangian case D.Talalaev considers: det{l — e~^''Tgi^_Yangian{z)) ([Ta04] formula 9 page 6). 



Remarks. 1) In |[GLZ03] ', '[Ko07B]', '[Ko07A]' the names "right quantum matrices"; in |[RRT02 
[[RRT05], .[LTOTj] the name "left quantum group" were used. We prefer to use the name "Manin 
matrices" . 



2) Our case is different from the more general one of |[GGRW02] where generic matrices with non- 
commutative entries are considered. In that case, there is no natural definition of the determinant 
(rather, it was proposed to build linear algebra in terms of elements A~^, elements of being 
called " quasi-determinants" , so that there are quasi-determinants instead of a single notion of 
determinant in the commutative case). Thus the analogues of linear algebra propositions are some- 
times quite different from the commutative case. Nevertheless some results of the above mentioned 
authors can be fruitfully applied to some questions here. 



See |[CT06J| [CM] 
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3) Those readers that are famihar with the i?-matrix approach to (quantum) integrable systems 
know that Lgi^_Gaudin{z) and Tgi^_Yangian{z) Satisfy quite different commutation relations (112II19P 
with " spectral parameter" , namely, in the first case we have linear i?-matrix commutation relations, 
while in the second one we have quadratic relations. It is quite surprising that Talalaev's introduction 
of dz converts both to the same class of Manin matrices. 

Moreover Manin's relations do not contain explicitly the "spectral parameter" z. Thus we can use 
simpler considerations (that is, without taking into account the dependence on spectral parameter) 
and apply our results also to the theory z-dependent Lax/transfer matrices. Another feature of inser- 
tion of dz is that the somewhat ad hoc shifts in the spectral parameter entering the formulas for the 
"qdet{T{z))" known in the literature now appear automatically from the usual column-determinant 
(see e.g. lemma O page [12]) . 

All the considerations below work for an arbitrary field of characteristic not equal to 2, but we 
restrict ourselves with C. 

Subsequent paper (s) [CFRuj will contain more proofs, background material, generalizations, e.g., 
to the (q)-super cases. 
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2 Manin matrices. Definitions and elementary properties 



We herewith recall definitions and results of '[Man87]', |Man88l IManQlj (with minor variations suited 
to the sequel of the paper) . We first remark that the notion (given in Definition [1]) of Manin matrix 
with elements in an associative ring TZ can be reformulated as the condition that 



Wp,q,k,l [Mpg,Mki] = [Mkq,Mpi] e.g. [Mn, M22] = [M21, M12], [Mu,M2k] = 



[M2i,M] 



ifcj 



Indeed, for g = / the above requirement yields that elements which belong to the same column of M 
commute among themselves. For g 7^ /, it is precisely the second condition of Definition [TJ 

Remark 1 Precisely these relations were explicitly written in |Man88j (see chapter 6.1, especially 
Formula 1). Implicitly they are contained in [Man87] - the last sentence on page 198 contains a 



definition of the algebra end{A) for an arbitrary quadratic Koszul algebra A. One can show (see the 
remarks on the page 199 top) that end{C[xi, Xn]) is the algebra generated by Mij. 



Remark 2 Relation with quantum group Funq{GL{n)) (see \Man88^ . ICFRuf ). Let us consider 
n = 2. The q-analog of the relations above is: M21M11 = gMiiM2i, M22M12 = gMi2M22, M11M22 = 
M22M11 + q~^M2iMi2 
Fung{GL{2)) (e.g. 



Man87 



gMi2M2i, this is precisely half of the relations defining the quantum group 
page 192, formula 3). This is true in general: q-Manin matrices 



are 
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obtained by half of the relations of the corresponding quantum group Funq{GL{n)). Conversely 
one can define Funq{GL{n)) by matrices M such that M and M* are simultaneously q-Manin 
matrices Man87] |Man88] . In present paper we consider q = 1 case. 

The proposition below gives a more intrinsic characterization of Manin's matrices (see |Man88j . 
|Man91j :) 

Proposition 1 Consider a rectangular n x m-matrix M, the polynomial algebra C[xi, ...,Xm] and 
the Grassman algebra CfT/^i, ...,'?/'„] T^-e. ipf = Ojipiipj = —ipjipi); suppose Xi and ipi commute with the 
matrix elements Mp^g. Consider the variables Xi, ipi defined by: 

- Mi,„, \ / \ / Mi,i ... Mi,„ \ 
... (^i,...,^J = (^i,...,^0 - , (3) 

- Mn,m / \Xm J \ M^,! ... M,,^ / 

that is the new variables are obtained via left action (in the polynomial case) and right action (in the 
Grassmann case) of M on the old ones. Then the following three conditions are equivalent: 

• The matrix M is a Manin matrix. 

• The variables Xi commute among themselves: [xi,Xj] = 0. 

• The variables ipi anticommute among themselves: ipiipj + = 0- 
Proof. It is a straightforward calculation. 

Let us present some examples of Manin matrices. 

Definition 2 A matrix A with the elements in a noncommutative ring is called a Cartier-Foata ( see 
ICF6S^ Fo79l ) matrix if elements from different rows commute with each other. 



Lemma 1 Any Cartier-Foata matrix is a Manin matrix. 

Proof. The characteristic conditions of Definition [1] are trivially satisfied in this case. □ 
The example below is related to our discussion of the Capelli identities (section I4.3.ip ). Let 
^ijiUij be commutative variables, and X, y be n x /c matrices with matrix elements Xij,yij. Also, 
let dx.dv he n X k matrices with matrix elements ir^, Tr^- Let z be a variable commuting with 
yij. The following 2n x 2k, and {n + k) x (n + k) matrices, with elements in the ring of differential 
operators in the variables X, Y are easily seen to be Manin matrices: 



2.1 The determinant 

Definition 3 Let M be a Manin matrix. Define the determinant of M by column expansion: 

rv 

cietM = def°'™M= ^(-1)" JJ M<,(i),„ (5) 

a£Sn i=l,...,n 

where Sn is the group of permutations of n letters, and the symbol r\ means that in the product 
Wi=i n ^o-{i),i one writes at first the elements from the first column, then from the second column 
and so on and so forth. 
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Lemma 2 The determinant of a Manin matrix does not depend on the order of the columns in the 
column expansion, i.e., 

crSSn i=l,...,n 

Proof. Since any permutation can be presented as a product of transpositions of neighbors (i, i + 1) 
it is enough to prove the proposition for such transpositions. But for them it follows from the equality 
of the commutators of cross elements (formula [1]) . □ 

Example 1 For the case n=2, we have 

det'"' « 5 V= «^ - c& " da - be. (7) 



^ c d 

Indeed, this is a restatement of the second condition of Definition [H 
2.1.1 Elementary properties 

The following properties are simple consequences of the definition of Manin matrix. 

1. Any matrix with commuting elements is a Manin matrix. 

2. Any submatrix of a Manin matrix is again a Manin matrix. 

3. If A, B are Manin matrices and Vi, j, /c, / : [Ajj, Bki\ = 0, then A + B is again a Manin matrix. 

4. If y4 is a Manin matrix, c is constant, then cA is a Manin matrix. 

5. If A is a Manin matrix, C is constant matrix, then CA and AC are Manin matrices and 
det{CA) = det{AC) = det{C)det{A). 

6. If A,B are Manin matrices and ^i,j,k.,l : [Aij.,Bki\ = 0, then AB is a Manin matrix and 
det{AB) = det{A)det{B). 

7. If A is a Manin matrix, then one can exchange the i-th and the j-th columns (rows); one can 
put i-th column(row) on j-th place (erasing j-th column(row)); one can add new column(row) 
to matrix A which is equal to one of the columns(rows) of the matrix A; one can add the i-th 
column(row) multiplied by any constant to the j-th column(row); in all cases the resulting 
matrix will be again a Manin matrix. 

8. If A and simultaneously A* are Manin matrices, then all elements Aij commute with each 
other. (A q-analog of this lemma says that if A and simultaneously A* are q-Manin , then A is 

quantum matrix: " RAA = AAR" (|[Man87J| |Man88j )). 



9. The exchange of two columns in a Manin matrix changes the sign of the determinant. If two 
columns or two rows in a Manin matrix M coincide, then det{M) = 0. 

Proof. The first assertion is obvious. For the second, (in the case of equal i-th and j-th 
column) use the column expansion of det{M), taking any permutation with first two elements 
i and j. One sees that det{M) is a sum of elements of the form {xy — yx){z), where x, y belong 
to the same column. Using the column commutativity [x, y] = we get the result. For the case 
of two coinciding rows the assertion follows at once, without the help of column commutativity. 
□ 

10. Since any submatrix of a Manin matrix is a Manin matrix one has a natural definition of minors 
and again one can choose an arbitrary order of columns (rows) to define their determinant. 
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2.1.2 The permanent 

The permanent of a matrix is a muhihnear function of its entries, similar to the determinant, with 
no (^—ly^dnier) j|-g (definition. We will make use of it in section [5l3] below. 

Definition 4 Let M be a Manin matrix. Let us define the permanent of M by row expansion^ : 



permM = pernf°'" M = ^ JJ M^,^©, (8) 

Lemma 3 The permanent of Manin matrices does not depend on the order of rows in the row 
expansion: ^ 

\/p e Sn perm''""M = T[ ^pW.-^bW) (9) 

crSSn i=l,...,n 

Proof. Since any permutation can be presented as a product of transpositions of neighbors {i,i + l) 
it is enough to prove the proposition for such transpositions. But for them it follows from the equality 
of the commutators of cross elements (formula [T]) . □ 

Example 2 / \ 

perm™- ( ^ ^ j arf + &c '™ da + cb, (10) 

Lemma 4 It is easy to see that the row-permanent of an arbitrary matrix M (even without conditions 
of commutativity) does not change under any permutation of columns. 



3 Lax matrices as Manin matrices 

In the modern theory of integrable systems, "Lax matrices" are also called "Lax operators" , "transfer 
matrices" , "monodromy matrices" for historical reasons. It is outside of the aim of the paper to make 
an attempt to provide a kind of formal definition of a Lax matrix for an integrable system. However, 
it is fair to say that there is a set of properties expected from a "good" Lax matrix; let us now recall 
those property relevant for our exposition. 

Let M be a symplectic manifold and Hi be a set of Poisson-commuting functions defining an 
integrable system on it (for simplicity, let us consider the Hamiltonian flow Xi of one of these 
functions, say Hi) A Lax matrix L{z) for (M, Hi) is a matrix, whose matrix elements are functions 
on M, possibly depending on a formal parameter z, satisfying the following two characteristic 
properties: 

1. The evolution of L{z) along Xi is of the form 

L{z) = [Mi{z),L{z)]. 



2. The characteristic polynomial det{X — L{z)) = "^ijHi jZ^X^ produces "all Liouville" integrals 
of motion, i.e. 

• Vz, j Hi j Poisson commute among themselves and with the given functions Hi 

Wk,i,j {Hk,HiJ = 0, yt,J,kJ {Hi,„Hk,i} = (11) 



^Remark the difference with the definition of the determinant, where one uses column expansion 
^"let us assume for simphcity of presentation that L{z) is just a polynomial function of the formal parameter z 
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• All Hi can be expressed via Hf^ i and vice versa. 



Poisson algebras of commutative functions on manifolds M are related to classical mechanics. As 
it is well known, in quantum mechanics one is lead to consider a family of non commutative, but 
associative algebras Fun[M)^] h is a. formal parameter in mathematics and Plank's constant in 

physics. For h = the algebra Fun{M)j^ coincide with the commutative algebra Fun{M) of functions 
on a manifold M, as well as Poisson brackets are related to commutators: [f,g] = h{f , g}'mod{h'^) , 



for f,ge Fun{M)^. (See, e.g., [K97]). 



The standard example is the algebra of functions on C^" - C[pi,gj] with the Poisson bracket 
{PiAj} = ^i,3i {Vi^Pj} = {liyQj} = whose quantization is the Heisenberg (Weyl) algebra generated 
by Pi, (ji and the relations [pi, qj] = h6ij, [pi,Pj] = [qi, qj] = 0. (We will usually put h = 1). 

Within this framework, on can look for a Quantum Lax matrix for a quantum systems as a matrix 
satisfying 

• L{z) is a matrix which matrix elements are elements from Fun{M)^, usually depending on a 
formal parameter z (thus, a quantum Lax matrix is a matrix with noncommutative elements). 



in the classical limit ^ ^ one has L(z) — > L(z 



\classical 



It is quite natural to look, in the integrable case, for a kind of determinantal formula: " det{\ — 
L{z))" = 'Ylk '^^Hk^z) to produce quantum integrals of motion: [Hk{z), Hiiu)] = 0, and, possibly, to 
satisfy other important properties (see ,[0X06]) . 



3.1 The Gaudin case 

The algebra of symmetries of the gln-Gandin ( |Ga76[ IGa83j ) integrable system is the Lie algebra 
gl[t]. The Lax matrix for the Gaudin system is a convenient way to combine generators of gl[t] 
(or its factor algebras) into one generating matrix- valued function (see formula [T8b below for gl[t] 
itself). Quantum commuting Hamiltonians (quantum conservation laws) arise from the maximal 
commutative Bethe subalgebra in U{gln[t])- Explicit efficient expressions for generators were first 
obtained by D. Talalaev. This subalgebra is actually an image of the center of Uc=crit{g^n[t,t~^],c) 
under the natural projection Uc=crit{gin[t,t~^],c) —>■ U{gln[t]) (see Appendix, section IHTTl) . This 
explains the mathematical meaning and the importance of the subject. 

Definition 5 Let R be an associative algebra over C. Let us call a matrix L{z) with elements in 
R{{z)) (i.e. L{z) e Main ® R® C{{z)))^ a Lax matrix of 5f/„-Gaudin type ijf: 

[Lij{u), Lki{v)] = —^{Lii{v)6jk - Lii{u)5jk - Lkj{v)5ii + Lkj{u)5ii). (12) 

More precisely this is "rational" (7/„-Gaudin type; as it is well-known,there are trigonometric and 
elliptic versions as well as generalizations to semisimple Lie (super)-algebras. 

We recall the matrix ("Leningrad's") notations: let 11 G Main ® Main be the permutation matrix: 
n(a ®h) = h®a, and consider L{z) L{u) e Main ® Mat„ ® i? (g) C((z)) (g) C((n)). Formula 

[T^ can be compactly written as follows: 

[L{z) 1,10 L{u)] = [-^, L{z) 01 + 10 L{u)], (13) 

z — u 

that is, a Lax matrix of 5fZ„-Gaudin type is a Lax matrix with linear r-matrix structure (the r-matrix 
being r = ^^). 

or L{z) G Main (S) R® €-{{z~^))] in many cases L{z) is just a rational function of z 
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Remark 3 FromD^) it follows that 



[Lij{u),Lki{u)] = -{duLu{u)5jk - duLkj{u)5ii). 
Here is our first main observation: 

Proposition 2 Consider dz ± L{z) G Main ® -R ® C((2;))[92], where L{z) is a Lax matrix of the 
gln-Gaudin type above; then: 



{dz — L{z)), {dz + L{z)Y are Manin matrices 
Proof. The proof is a straightforward computation. 
Example 3 Let us show this property in the example of 2 x 2 case: 
dz-Lii{z) -Li2{z) 



(14) 



a - Liz) 



-L2i{z) dz-L22{z) 



Column 1 commutativity : [dz — Lii{z), —L2i{z) 



-L' 



21^ 



[L^^{z),L2i{z)]=Q, 



(15) 



Cross term relation : [dz - Lii{z), dz - L22{z)] = -L'22{z) + Ll^^{z) = [L2i{z), Li2{z)]. 

The following well-known fact easily follows from [T31 

Proposition 3 Let L{z) and L{z) be Lax matrices of the Gaudin type with pairwise commuting 
elements (i.e. \/i, j, k,l, z,u : [Lij{z), Lki{u)] = 0) then L{z) + L{z) is again a Lax matrices of the 
Gaudin type. 



Remark 4 The classical counterpart of the commutation relations (fT2|) - namely, with Poisson 
brackets replacing commutators) and of the corresponding Lax matrices are associated with a huge 
number of integrable systems. In particular, among integrable systems leading to Lax representations 
of Gaudin type we can mention: the Neumann's geodesic motion on ellipsoids and the Neumann- 
Rosochatius systems, the ^f/^-Manakov-Mishenko-Fomenko tops (see, e.g., [[AHH96] ) to the so(3) 
case); the Lagrange top and the Clebsh system (see, e.g., [PR07]); "polynomial matrices" (and so 



finite gap solutions of KdV-type equations, although non-trivial reductions of the r-matrix bracket 
have to be considered) |Mum84] . [RS79], [AvM80], [Be90]i the general systems studied in [AHP887i 
the Jaynes-Cummings-Gaudin system, jBTa07] • as well as Hitchin's system on singular rational curves 
[[CT03B] . Bending fiows, Gelfand-Ceitlin, Jucys-Murphy integrable systems are degenerations of the 
Gaudinlystem (I [FM03A] | [FM03B] | [CFRy07] I) . 



Moreover, Schlesinger system of isomonodromy deformations and Knizhnik-Zamolodchikov equa- 
tion are non-autonomous versions of the classical and, respectively, quantum Gaudin Systems [DM03] , 
[[Re92], It is well-known that Nahm's monopole equations can be rewritten as a Gaudin system (see 
e.g. |[GG M97] (section 6.1 page 24), where applications to Seiberg-Witten theory of N=2 supersym- 
metric Yang-Mills are also discussed). In the infiuential paper FFR94 (see, also the survey Fr95 ) 
the relation with the geometric Langlands correspondence was discovered. It might also be noticed 
that a great number of systems in condensed matter theory are related to Gaudin-type systems, see, 
e.g. 



OSDR04 



In this paper, we shall concentrate on the following few examples of Gaudin type Lax matrices. 
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The simplest example Let K be an arbitrary constant matrix, and n, /c G N, and zi^...^Zk 
arbitrary points in the complex plane. Consider 



i—l,...,k y Quyi / 

where pij, g^j-, z = 1, n; j = 1, k are the standard generators of the standard Heisenberg algebra 
\Pi,jAk,i\ = 5i,k^j,h [Vij^Vk^i] = Vli,jAk,i\ = 0. Also, Q,P are n x fc-rectangular matrices with elements 

One can see that this Lax matrix satisfies relations [13] f[T2l) . and so by proposition [2] above 
[dz — L{z)), [dz + L{z)y are Manin matrices. 

The standard example Consider gin ® ■■■®gln and denote by e\i the standard basis element from 
the z-th copy of the direct sum gl^ © ... © gin- The standard Lax matrix for the Gaudin system is: 

Lgl„—Gaudin standardi^Z^ ^ ^ 1 " I ("^''') 

1=1,. ..,k 

Lgl^-Gaudin standard{z) E Matn ® U {gin © ... © gin) ® '^{z) ■ 

One can see that this Lax matrix satisfies relations [T3] (1121) . and so by proposition [2] above 

dz — Lgi^^-Gaudin standard{z), {dz + Lgi^^^^Gaudin standard{z)y are Manin matricCS. 

The gin and gln[t] examples. Consider the Lie algebra gin, with e^j its standard linear basis; 
consider the polynomial Lie algebra gln[t]/t^, k = l,...,oo. The Lax matrices for gin and gln[t]/t^ 
are the following: 

... ei^n \ 1 / ei,if-^ ... ei,„f-^ 



ei,i - 


p* 




^n,n 





p f 



i-1 



Lgi^{z) e Matn®U{gln)®C{z)] Lgi^^t]{z) e Matn^U{gln[t])^C{{z-'^)), where U{g) is the universal 
enveloping algebra of 5^. Note: zLgi^{z) coincide with the expression (12) section 2.1 [KirOO]. 

One can see that these Lax matrices satisfy relations [13] f[T2l) . and so by proposition [2] above 
dz - Lgi^{z), dz - Lgi^^t]{z), {dz + Lgi^{z))\ {dz + Lgi^^t]{z)y are Manin matrices. 



3.2 The Yangian (Heisenberg chain) case. 

The Algebra of symmetries of the Heisenberg XXX spin chain, the quantum Toda system and several 
other integrable systems is a Hopf algebra called Yangian. It was implicitly defined in the works of 
Faddeev's school; a concise mathematical treatment and deep results were given in jDr85j (see also 
the excellent more recent surveys [Mo02]^jMNO94] ). It is a deformation of the universal enveloping 
algebra of gl\t]. A Lax matrix of Yangian type is a convenient way to combine generators of the 
Yangian (or its factor algebras) into one generating matrix-valued function. Following standard 
notations we will write T{z) ("Transfer matrix") instead of L{z) in the case of the Yangian type Lax 
matrices. 
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Definition 6 Let TZ he an associative algebra over C Let us call a matrix T{z) with elements in 
TZ{{z)) (i.e. T{z) G Main ® 7^ ® C{{z)))^ a Lax matrix of Yangian type iff: 

[Tij{u),Tki{v)] = -^{Tkj{u)Tu{v) - nj{v)Tu{u)) (19) 

(See |[Mo02]' page 4 formula 2.3. More precisely this is the case of the Yangian, there are 
generalizations to the semisimple Lie (twisted)- (super )-algebras). 

In matrix (Leningrad's) notations we have, with 11 G Mat„ ® Mat„, the permutation matrix: 
n(a®6) = 6®a, andT(2)®l,l®T(M) G Ma^„®Ma^„®7^®C((^))®C((M)), R{z-u ) = (1® 1-^) 
that formula [T^] can be written as follows, as a quadratic i?-matrix relation: (see [Mo02 J page 6 
proposition 2.3 formula 2.14). 

(1 ® 1 - -^)(T{z) ® 1) (1 ® T{u)) = (1 ® T{u)) {T{z) ® 1)(1 ® 1 - (20) 
z — u z — u 

Or shortly: R{z - u)T{z)T{u) = T{u)T{z)R{z - u). (21) 

Our second main observation is: 

Proposition 4 IfT{z) is a Lax matrix of the Yangian type then 

T(z)e~^'', {e^''T{z)Y are Manin matrices. (22) 

Here e^^'^T^z) G Main ® Y{gQ ® C[[l/2, e"'^^]], (see definitionlE above). Actually Vc G C : r(2 + c) 
is a/so of the Yangian type, so T{z + c)e-^'~ = e-'^^T(z + c + 1), (e^'-T(2 + c))* = {T{z + c + l)e'^^)* 
are Manin matrices. Efl 

Proof. As in the Gaudin-type case, it follows from a straightforward computation. Actually it is 
known that any k x k submatrix of the Yangian type matrix is again a Yangian type matrix of 
k X /c-size, so essentially it is necessary to check only 2x2 case. 
The following well-known fact easily follows from| 



Proposition 5 Let T{z) and T{z) be Lax matrices of Yangian type with pairwise commuting ele- 
ments (i.e. "^ijjjkJjZjU : [Tij{z),Tki{u)] = 0); then the product T{z)T{z) is again a Lax matrix of 
Yangian type. 



Lemma 5 Observe that det™'"'"''(T(z)e-^^) = grfet(T(z))e-"'='% where 

qdet{M{z)) = Eaes„ Hfc ^^(fc),fc(^ - k + 1) (see e.g. formula 2.24 page 10 |[Mo02]D . From 2.25 
loc.cit. one sees: det^"^{e^^T{z)) = e'^^'^ qdet{T (z)) . So "shifts" in z in formulas for qdet now appear 
automatically from the usual column/row determinants. 



Let us herewith list a couple of remarkable examples of Yangian-type Lax matrices. 

^'^ or T{z) G Main ® 7?. ® C{{z~^)); in many cases T{z) is just a rational function of z 

1 2 

^^Observe an analogy. It is known ( [Man88j ) that, if M is matrix defining Funq{GLn) i.e. it satisfies RM Ad = 

2 1 

MM R without spectral parameter, then Ad and simultaneously Af* are q-Manin matrices. Here we have: T{z) 

12 2 1 

satisfies R{z — u)T{z)T{u) = T{u)T{z)R{z — u) with spectral parameter and proposition claims that e~ "T^z) and 
(e^^r(z))* are simultaneously Manin matrices. Both cases are related to GL„, extension to the other Lie groups is 
not known for the moment. 
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The Toda system Consider the Heisenberg algebra generated by ^ = 1, and relations 

= Sij, [pi,Pj] = [qi,qj] = 0. Define 



TtoM = n ( '7 ) (23) 



-e''' 

i=l,...,n 

One can see that this Lax matrix satisfies relations [2U] f[T9]) . and so by propositionlHabove e~^^TToda{z) 
is a Manin matrix. One can easily see that, with the identification qn+i = qi, the coefficient C„_i of 



z 



n-2 



oiTrTToda{z) equals Y.k<i<nPkPi-Yji=i,...,n^'^" ''"^^^ and the coefficient Cn-i of ^ is - YH=iPi- 
Thus |C^_i — Cn~2 is the physical Hamiltonian (i.e., the energy) of the periodic Toda chain. 

The Heisenberg's XXX system Consider the Heisenberg algebra generated by [Pi,jiqk,i\ = 
Si,kSj,i, [Pi,j,Pk,i] = [qij, qk,i] = 0; the quantum Lax matrix for the simplest case of the ^fZ^- Heisenberg 
spin chain integrable systems can be given as follows: 



-XXX "simplest" 




(24) 



One can see that this Lax matrix satisfies relations [20] f|T9l) . and so by proposition H] above 
e~^''Txxx{z) is a Manin matrix. 

Consider gin © ... ®gln and denote by e^i the standard basis element of the i-th copy of the direct 
sum gin © ... © gin- The standard Lax matrix for the GL(A^)-Heisenberg's XXX system is: 

/ 1 / ell ••• eln \ \ 

Tgl^-XXX-standard{z) = Y\_ I ^nxn H [ j j (25) 



Tgl„-XXX-standard{z) G Motn ® U {gin © ... © gL) ® C{z) 

One can see that this Lax matrix satisfies the relations [20] ([T^ . and so by proposition H] above 



e ^''Tgi^_xxx-standard{z) is a Maulu matrix. 



4 Algebraic properties of Manin matrices and their appli- 
cations. 

In this section we will derive a few less elementary properties of Manin matrix, and give applications 
thereof to integrable systems. 

4.1 Cramer's formula 

Proposition 6 \Man88^ Let M be a Manin matrix and denote by M"''^^ the adjoint matrix defined in 
the standard way, (i.e. M^f = (— l)^"'"'(iet'^°'"'""(M;fc) where Mik is the (n — 1) x (n — 1) submatrix of 
M obtained removing the l-th row and the k-th column. Then the same formula as in the commutative 
case holds true, that is, 

W'^m = def'''""*"(M) Id (26) 

//M* is a Manin matrix, then M°"^^ is defined by row-determinants and MM"-^^ = (iet™'"'""'(M*) Id = 
cZet™"'(M) Id. 
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Example 4 In the 2x2 case we have: 

f d —h\fah\_fda — hc dh — hd \ _ f ad — ch 
\ —c a J \ c d J y —ca + ac —cb + ad J y ad — cb 

where the characteristic commutation relations of a Manin matrix have been taken into account. 

Proof. One can see that the equality: V i : {M"-'^^ M)i^i = det^°\M), follows from the fact 
that det'^°\M) does not depend on the order of the column expansion of the determinant. This 
independence was proved above (lemma [2]). Let us introduce a matrix M as follows. Take the matrix 
M and set the i-th column equals to the j-th column; denote the resulting matrix by M . Note that 
det^°\M) = precisely gives {M°''^^ M)ij = for i ^ j. To prove that det'^°\M) = we argue as 
follows. Clearly M is a Manin matrix. Lemma [2] allows to calculate the determinant taking the 
elements first from i-th column, then j-th, then other elements from the other columns. Now it is 
quite clear that det'^°\M) = 0, due to it is the sum of the elements of the form [xy — yx){z) = 0, 
where x, y are the elements from the i-th and j'-th of M, so from j'-th column of M. By Manin's 
property elements from the same column of a matrix M commute, so xy — yx = 0, so det'^°\M) = 0. 
□ 

Remark 5 The only difference with the commutative case is that, in the equality (!26|) the order 
of the products of M"-'^^ and M has to be kept in mind. 



4.1.1 Application to the Knizhnik-Zamolodchikov equation 

We can give a very simple proof of the formula relating the solutions of KZ with coupling constant 
K = 1 to the solutions of the equation defined by Talalaev's formula: 

det{d, - L{z))Q{z) = 0. 

This result was first obtained in '[CT04]' in a more complicated manner (see alsol[CT06b]i). 

The standard gf/^-KZ-equation (actually a system of equations) jKZ84]j for an unknown function 
^(^0, zi, Zk) is the following. Let Vi,i = 0, k be representations of gin, "^{zq, zi, Zk) is 
Vo ® Vi ® ... ® Vfc valued function of Zp, p = 0, k. Let Cat be the standard basis in gin, the 
corresponding basis in gin © ... © gin (upper index i is a number of copy of gin)- Let vr denote the 
representation of gin © ... © gin in Vq © Vi ® ... ® V^. 

KZ system. j-th equation: \ - k V ^Lj±!1^^2k>^l!d^ ^(^) = q. (28) 

j = 0, k, and k G C is some constant. 

Let us restrict to the case Vq is and denote zq by z. Let us now consider only one equation, 
say the first and observe the following relation with the Gaudin Lax matrix: 



=i...fc 




So we came to: 

Observation. ('[CT06b]l) The equation {dz — 7r{LGaudin{z))'^ (z) = is precisely the particular case 
of the KZ-equation0 



^^From the point of view of the Gaudin system the first space Vq = C" - is the " auxihary space" ; and the tensor 
product of the other spaces Vi (Si ... ^ Vk is the "quantum space" (using integrabihty theory language). The auxihary 
and quantum spaces play the very different roles in integrabihty. But in KZ-equation all the spaces 1^, i = 0, ... play 
the same role. This uncovers unexpected symmetry between auxiliary and quantum spaces. It is different from the 
other relations between KZ-cquation and the Gaudin model (e.g. FFR94 J. 
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Here Lcaudiniz) = Xli=i k '' is the standard Lax matrix of the quantum Gaudin system 

(see formula [T7] page [T2D considered in a representation vr; {n,Vi ® ... Cg) V^) is a representation of 
{gin © ... © gin)', ^(^) is a C" © Vi ® ... © V^- valued function, so that its components "^{{z) are 
Vi © ... © Vfc-valued functions; Eat G Main and Cat G C U{gln) are the standard matrix units. 

Proposition 7 Lei 6e a solution of the KZ-equation with k = 1 [29\] . then: 

yt = 1, n Tl{det{d, - LGaudin{z)))^i{z) = (30) 

Proof. The adjoint matrix {dz — Lcaudiniz))"''^^ exists as we discussed above (here we use n = 1), 
so: 

TT{dz - LGaudiniz))<il{z) = 0, ^ 7r((9, - ^Ga^dm (2^) )'"*>(9, - LGaudiniz))^ (z) = 0, 

hence 7r{det{dz - LGaudin{z)))Id "^{z) = 0, (31) 
hence Vz = 1, n 7i{det{dz - LGaudiniz)))"^ i{z) = 0. 

□ 



The equation det{dz — L{z))Q{z) = should be seen as generalized Baxter's T — Q- 
CT06] ), it is well-known that Baxter's equation plays a crucial role in the solution of 



Remark 6 

equation (see 

quantum systems. In particular it is important to establish that its solutions are rational functions 
(see MTV05 ). The proposition above relates this question to the rationality of KZ-solutions which 
is a more natural question (see |[Sa06]| ). 



4.2 Inversion of Manin matrices 

Theorem 1 Let M be a Manin matrix, and assume that a two sided inverse matrix exists 
(i.e. M~^M = MM~^ = Then is again a Manin matrix. 



Proof. The proof of the assertion consists in a small extension and rephrasing of the arguments 



in the proof of Lemma 1 page 5 of [BT02 



We consider the Grassman algebra C[ipi_, V'n] (i-e. ipf = Oyipi'ipj = —ipj'ipi), with ipi commuting 
with Mp g. Let us introduce new variables ipi by: 

Mi,i ... Mi,„ 

(^i,...,7/.„) = (V'i,...,V'n,) I ... |. (32) 

Mn,l ... Mn,n 

It is easy to see that the Manin relations [M^-^Mj^j^] = [M^^^, M^-^] we have to prove follow from 
the equahty 

iW,-^M-Vi A ... A ^„ = M-^M-^ij^ A ... A^n + 

_ _ i—th place k—th place _ 

+ (rfet(M))"Vi A^2 A ...A ipj A.. .A ipi A...Aipn- (33) 

(Here assume i k; otherwise the desired result is a tautology). Let us prove (!33|) . By the definition 
of the ipi^s, {ipiy ■■■,4'm) = ("^1, ■■■,4'n)M, multiplying this relation by on the right we get 

{i>i,...,i>n)M-^ = {^P,,...,^|Jn), ^ ^^,M;i^=i)i. 
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i-th place k~th place 

We multiply this relation with the n — 1-vector ipi Aip2 ^ •••A ipj A... A empty A... A ipn, and 
use the fact that Vm : ip'^ = 0, to get 

„ _ „ _ i—th place k-th place 

{ipiM-^ + ijkM^^^ -iji)^iAij2A--A i^j A. ..A empty A...AV^„ = 0. 

_ j—th place _ 

Now, using: det{M)M~^ tpi A ... A ipn = i'l /\ •••A tpi A... A tpn we get the desired result. □ 



Remark 7 The paper [Ko07B] (see section 5 page 11 proposition 5.1) contains a somewhat 



weaker proposition for matrices of the form 1 — tM, where t is a formal parameter. It states that for 
(l-tM)-^ the cross-commutation relations holds true. The "column commutation" property is not 
considered there. The equality (iet™'*""'^((l — tM)~^) = ((iet™'"'"'^(l — tM))~^ is contained in theorem 



5.2 page 13 of |[Ko07B [ The proofs presented there are based on deep combinatorial arguments. 



4.2.1 Application to the Enriquez - Rubtsov - Babelon - Talon theorem. 



Let us show that the remarkable theorem '[EROl]' (theorem 1.1 page 2), '[BT02]' (theorem 2 page 4) 
about "quantization" of separation relations - follows as a particular case from theorem [1] above. 

Let us briefly recall the constructions of this Theorem (following, for simplicity [BT02]). Let 
{<yi, (3i}i=i,...,g be a set of quantum "separated" variables, i.e. satisfying the commutation relations 

[ai,aj] = 0, [/3i,pj] = 0, [ai, (3j] = f{ai, f3i)6ij, i,j,= l...,g. 

satisfying a set of equations (i.e., quantum Jacobi separation relations) of the form 

a 



^ Bj{ai, pi)Hj + Bo{ai, pi) = 0, i 



.9, 



(34) 



for a suitable set of quantum Hamiltonians Hi, . . . , Hn. One assumes that some ordering in the 
expressions Ba,a = 0,. . .g between has been chosen, and that the operators Hi are, as it is 

written above, on the right of the Rj. Also, the Ba are some functions (say, polynomials) of their 
arguments with C-number coefficients. Then the statement is that the quantum operators Hi, ... , Hn 
fulfilling ( IMI) commute among themselves). 

Our proof starts from the fact that one can see that equation fl3^ can be compactly written, in 
matrix form, as 

B-H = ~V, with Bi, = Rj{ai, (3i), Vi = Boiai, A), 



and thus (making contact with the formalism of [EROl] ) one is lead to consider the gx (g + 1) matrix 



A 




Thanks to the functional form of the -Bjj's and of the V^'s, this matrix is a Cartier-Foata matrix (i.e., 
elements form different rows commute among each other), and hence, a fortiori a Manin matrix. 

Given such a. g x {g + 1) Cartier-Foata matrix A, we consider the following {g + 1) x [g + 1) 
Cartier-Foata (and hence, Manin) matrix: 



A 



( 1 

Vi Bi,i 



\ 



9,1 



B. 



9,9 



J 
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Now it is obvious that the solutions Hi of equation (15^ are the elements of the first column of the 
inverse of A, and namely 

Hi = = 1, . . . ,5f. 

Since A is Cartier-Foata, its inverse is Manin, and thus the commutation of the His can be obtained 
from theorem [TJ □. 



Remark 8 For the sake of simplicity, we considered, as in [BT02] , an arbitrary system for which 
a quantized spectral curve exists, and the commuting Hamiltonians can be written by the formula 
above in terms of the so-called separated variables (see section [6]). However, it should be noticed 
that the theorem hold in a more general (in classical terms, Stackel or Jacobi) setting, where the 
functional dependence of the matrix elements Bij{ai, Pi) and Vi{ai, (3i) on their argument may depend 
on the index i. 

4.3 Block matrices and Schur's complement 

Theorem 2 Consider a Manin matrix M of size n, and denote its block as follows: 

= { r n 

\ ^n—kxk -'-^n—kxn—k / 

Assume that M,A,D are invertible i.e. 3M~^,A~^,D^^ : A'^A = AA'^ = 1, D-^D = DD'^ = 1, 
M~^M = MM^^ = 1. Then the same formulas as in the commutative case hold, namely: 

def''\M) = def°\A)def°\D - CA^^B) = det'^\D)def''\A - BD-^C), (36) 

and, moreover, the Schur's complements: D — CA^^B, and A — BD^^C are Manin matrices. 

Sketch of the Proof. The full proof of this fact will be given in [CFRu] . Actually it is not 
difficult to see that A — BD~^C, D — CA~^B are Manin matrices from theorem [1], so basically we 
need to prove the determinantal formula fl5^ . To this end, one proves the 

Lemma 6 Let M he a Manin matrix of size n xn. Let X be a k x [n — k) matrix, for some k, with 
arbitrary matrix elementJ^. Then: 

^^^colun^n^^ f Uxk ^^xn-k \ ^ ^^^column ^37) 
\ 'Jn—kxk ^n—kxn—k J 

Now, the desired result follows now from the equality: 

A B\fl -A-^B\^fA \ 

C D J \ 1 J \C D- CA-^B J ■ ^ > 

Remark 9 1) Alternatively, one can prove fl5^ using the Jacobi type theorem by I. Gelfand, V. 
Retakh (see e.g. [GR97 theorem 1.3.3, page 8) with similar arguments to those used by D. Krob, B. 



Leclerc KL94 (theorem 3.2 page 17). 

2) The Jacobi ratio theorem was proved for Manin matrices of the form 1 — tM, t is a formal 
parameter, in the remarkable paper [Ko07B](see theorem 5.2 page 13). It is actually equivalent 
to our formula [3S] for matrices of such a form. As in the case of the matrix inversion formula, 
column commutativity is not considered, and the proofs therein contained are based on combinatorial 
properties. 



^No conditions of commutativity at all are required. 
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4.3.1 Application. MTV-Capelli identity for gl[t] (for the Gaudin system) 

In l[MTV06]' a remarkable generalization of the Capelli identity has been found. Using the beautiful 

insights contained therein, we remark that it follows from theorem [2l 

Consider C[pij, qij], i = 1, n;j = 1, k, endowed with the standard Poisson bracket: 
{Pi,jyQk,i} = Si,kSj,h {Pi,j^Pk,i} = {^ij! = 0. Cousider its quantization i.e. the standard 

Heisenberg (Weyl) associative algebra generated by Pij, qij, i = 1, ...,n;j = 1, ...,k, with the relations 

\Pi,j,qk,i] =Si,kSj,i, [Pi,j,Pk,i] = [qi,j,qk,i] = O. Define n xk-rectangulai matrices Q classical, Pciassicai,Q,P 
as follows with: Qij = qij, Qij = qij, Pij = pij, Pij = pij. Let Ki,K2 he n x n, k x k matrices 
with elements in C. 

Let us introduce the following notations: 

^.uantum^^^ = i^l + Q{Z - K,)-'P\ L^'™'(Z) =K, + Qclass^cal{z - K^T' P^ssical (39) 

These are Lax matrices of Gaudin type i.e. they satisfy the commutation relation f[T2]) . 
Let us reformulate the result of MTV06 in the following simple form: 

Proposition 8 Capelli identity for the gl\t\ (Gaudin's system) 

Wick{det{X - L^i'^^^'^^^lz))) = det{d, - ^^"""^""^(z)) (40) 
Here we denote by Wick the linear map: C[\,Pij,qij]{z) C[dz,Pij,qij]{z), defined as: 

Wzck{f{z))^l[q^^^l[t^) = f{z)d:llq;^l[f^ (^1) 

ij ij ij ij 

i.e. that, from a commutative monomial, makes a noncommutative monomial according to the rule 
that any q placed on the left of any p, and the same for z and dz- It is sometimes called "Wick or 
normal ordering" in physics. 

Proof. Following MTV06]^ we consider the following block matrix: 

MTV = ( ' " ^2 ^ ^* 1 (42) 

It is easy to see that MTV is a Manin matrix(indeed it is strictly related to the Manin matrices 
briefiy introduced in (jlj)) ; also, as it was observed by Mukhin, Tarasov and Varchenko, the Lax 
matrix of the form 0391) appears as the Schur's complement "D — CA~^B" of the matrix MTV: 



z 



By theorem [2] we get that 

det'oiumn^^rpy^ = det{z - K2)rfet™'""^"(9, - K2 - Q{z - lUY^P'). (43) 



This a form of the Capelli identity as explained in [MTV06]'. 

In order to arrive at formula fHOl) above, we just remark that in det'^"''^"^"' (MTV) all variables z, qij 
stands on the left of the variables dz,Pij- This is due to the column expansion of the determinant, 
where at first appears the first column, then the second, and so on an d so forth. Actually, the 
operators z, q^j stand in the first ra-th columns of MTV matrix, while the operators dz,Pij are in the 
m rightmost columns. So dividing fH51) by det{z — K2) we obtain fl4UI) . DSee section 18.3.21 and the 
references quoted therein for the relation with the classical Capelli identity. 
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5 Spectral properties of Manin matrices and applications 
5.1 The Cayley- Hamilton theorem 

The Cayley-Hamilton theorem, (i.e, that any ordinary matrix satisfies its characteristic polynomial) 
can be considered as one of the basic results in linear algebra. The same holds - with a suitable 
proviso in mind - for Manin matrices. 

Theorem 3 Let M be a n x n Manin matrix. Consider its characteristic polynomial and the coeffi- 
cients hi of its expansion in powers oft: det'^"'-^"^^ [t — M) = Yli=o n^«^*'' ^hen 

J2 h^M' = I.e. rfef°'""^"(t-M)|;it? = 0. (44) 

If M* is a Manin matrix, then one should use left substitution and row determinant: 
def°'"{t - M)|fJ{/"''^***"*'= = 0. 

Proof. It was proved in Proposition [6] that there exists an adjoint matrix (M — t Jd)"'^^ , such 
that 

(M - t Idf^^iM -t Id)= det{M - t Id)Id. (45) 

The standard idea of proof is very simple: we want to substitute M at the place of t; the LHS of this 
equality vanishes manifestly, hence we obtain the desired equality det{M — t Id)\t=M = 0. The only 
issue we need to clarify is how to substitute M into the equation and why the substitution preserves 
the equality. 

Let us denote by Adjk{M) the matrices defined by: J2k=o,...,n-i^djkiM)t^ = {M -t Idf^K The 
equality above is an equality of polynomials in the variable t: 

(^ Ad3k{M)t''){M -tId) = J2 AdUM)Mt' - ^ AdUM)t''+' = (46) 

= det{M -t Id)=J2 hkt''. (47) 

k 

This means that the coefficients of t' of both sides of the relation coincide. Hence we can substitute 
t = M in the equality, substituting "from the right": 

Y.AdjkMM'^ -Y,^d3kM^+^ = Y.hkt\^^^. (48) 
The left hand side is manifestly zero, so we obtain the desired equality: Ylik=a n^kM^ = 0. □. 



Remark 10 See Mo02 1 (section 4.2), GIQPS05], and references therein for other CH-theorems. 



5.1.1 Example. The Cayley-Hamilton theorem for the Yangian 

Let us consider the Lax matrix Tgi^_Yangian{z) (or T{z) for brevity) of Yangian type (see section! 
The matrix T(z)e~^^ is a Manin matrix. Let us derive a Cayley-Hamilton identity for T{z) from the 
one for the Manin matrix T{z)e~^''. 

Definition 7 Let us define the "quantum powers" T^^{z) for the Yangian type Lax matrix as follows: 
T^]{z) =^ T{z + p) -Tiz+p-l) ■ ...-Tiz + l) = eP^^{T{z)e-^^y (49) 
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Denote by QHi(z) coefficients of the expansion of the "ordered" characteristic polynomial in 
powers of e^"": 

cief'''"™"(t-T(2)e-^^) = t^QHk{z)e^^-''^^% (50) 

k=0,...,n 

Note that the operator dz does not enter into the expressions QHi{z). Explicit formulas for QHi{z) 
are obviously the following: 

QHr,{z) = 1, Qi/„_i(z) = -Tr{T{z)), QHo{z) = (-1)" qdet{T{z)), (51) 

gi^„_,(2) = (-i)(«-) En<...<n qdet{T{z),,,....,:). (52) 

In words, the QHi{z)^s are the sums of principal q- minors of size i (in complete analogy with the 
commutative case, modulo substitution of minors by q-minors). Recall that qdet{M{z)) is defined 



by the formula qdet{M{z)) = Xl<Te5„ Ilk Ma{k),k{^ - k + 1) (see e.g. formula 2.24 page 10 |[Mo02][ ). 



Proposition 9 The following is an analogue of the Cayley- Hamilton theorem for Yangian type Lax 
matrices: 

J2 QHk{z + n)T^''\z) = 0. (53) 

k=0,...,n 

Proof. The proof is very simple: one uses Cayley-Hamilton theorem [3] for the Manin ma- 
trix T{z)e~^% and than excludes powers of e^^ from the identity!^. We can proceed as follows. 
By theorem [3] it is true: 

h,{z,dz){T{z)e"'^)' = 0, Y h,{z,dz)e"''^T^'\z) = 0, (54) 

k=0,...,n k=0,...,n 

where hk{z,dz) are defined as det'^°^^^'^{t — T{z)e~^'') = J2i=o n^i(^^^z)^^- definition of 

QHi{z) (formula USD one has: hi(z,dz) = QHi{z)e^^~^^^' , substituting this into EH one obtains: 



0= ^ hk{z,dz)e-''^^T^''\z) = J2 QHk{z)e^^-'''^^'e-^^^T^^\z)= (55) 

k=f),...,n k=0,...,n 

= J2 QHk{z)e-''^^T^'\z) = e-""^' Yl QHk{z + n)T^''\z), (56) 

fc=0,...,n k=0,...,n 

So we conclude that T.k=o,...,nQHk{z + n)T^^\z) = 0. □ 

Remark 11 The quantities QH^^z) were defined in |KS81] (formula 5.6 page 114) by means of 
different arguments. The idea of obtaining them as coefficients of a suitable characteristic polynomial 



is due to D.Talalaev Ta04 



5.2 The Newton identities and applications 

As it is well known, the Newton identities are identities between power sums symmetric functions and 
elementary symmetric functions. They can be rephrased as relations between TrM^ and coefficients 
of det{t + M) for matrices with commutative entries. We will show that the same identities hold true 
for Manin matrices and present applications. 



^^we slightly changed the definition of QHi(z) comparing with the version 1, and correct the misprint in formula lSOl 
^^One can easily exclude dz in the Yangian case, but it is not clear at the moment how to do it in the Gaudin case. 



Nevertheless, a similar identity can be proved in for the Gaudin case also: [CT06b] (by different methods) 
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Theorem 4 Let M be a Manin matrix, and denote by the coefficients of the expansion in t of 
the determinant oft — M, i.e., rfer°'"""^(t — M) = hkt^ ; conventionally, let hk = for k < 0. Then 
the following identity holds: 

afCiet™'"'"'^(t-M) = ^(c/ef°'"'""(t-M)) ^ Tr{M/tf, (57) 

/c=0,..-,oo 

<^ \/k : -oo < k<n it holds: khk = ^ hk+iTr{M)\ (58) 

i:max{0 ,—k)<.i<n—k 

where n is the size of the matrix M . If is a Manin matrix, then dtdet^°^{t — M) = {det'^°^{t — 
M)) {lltT.k=o,...,ooTr{M/tf). 

So these identities are identical to those of the commutative case, provided one pays attention to 
the order of terms: hiTrM^ if M is a Manin matrix {TrM^hi if M* is a Manin matrix). Obviously 
enough, this difference is due to the absence of the commutativity. 

Sketch of the Proof. First observe the following simple property, whose proof is straightfor- 



ward^*^'- 



Tr{t + Mf^^ = dtdet{t + M). 

So we consider 



l/t Tr{{-M/tf) = Trj^ = Tr({def°^'''''''{t + M))-^{t + Mf'^^'^ 

fc=0,...,oo 

= (c^ef °'""^'^(t + M))-^Tr{t + Mf'^^ = (rfef "'"'""(t + M))-^dtdef°^'''^''{t + M). (59) 



Substituting —M instead of M one obtains the result. The case where M* is a Manin matrix is 
similar. □. 

Remark 12 While working at the proof of this proposition we were informed by P. Pyatov (un- 
published, based on Gurevich, Isaev, Ogievetsky, Pyatov, Saponov papers |GIOPS95] - |GIOPS05j ) 
that the Newton identities holds true as well as their q-analogs. Our proof is different and very 
simple, but it is a challenge how to generalize such an argument to the case g 7^ 1. 

5.2.1 Quantum powers, ZUcrit{gl'n\tit~^]^c), and new quantum Gaudin integrals 



We recall from [CT06b] a possible definition of "quantum powers" of the Lax matrix of the Gaudin 
model and apply the quantum Newton relations above to prove the commutativity of their traces. 
Hence we obtain new explicit generators in the center of Uc=crit{g!'n[t, t~^], c) and in the commutative 
Bethe subalgebra in U{gln[t])) (see Appendix (section [HH]) for a reminder); in a physicist's language 
we give explicit formulas for quantum conservation laws for the quantum Gaudin system. They are 
quantum analogs of T^i^i^dassicai^^^Y ^ and the result below exhibits the appropriate corrections for 
traces of powers which preserves commutativity. 

Definition 8 The quantum powers of Gaudin type Lax matrices (definition\^ formula lT^) are defined 
inductively as follows: 

lM(z) = Id, L[^1(z) = L^'-^\z)L{z) + {L^'-^\z)y (60) 



Actually, the equality holds for any matrix M , provided one consistently defines the determinant and the adjoint 
matrix. 
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Here (L^^~^^{z)y is the derivative with respect to z of (L'*"^] (2;)). ^^(z) are noncommutative ana- 
logues of the Faa di Bruno polynomials ( |Di03] section 6A, page 111). (Remark that in the commu- 
tative case L{z) and L{z)' commute). The binomial type formula below is due to D. Talalaev: 

(9. - L{z)y = J2 (-1)" ( ' ) dt'^L^^\z). (61) 

p=0,...,i V ^ / 



Theorem 5 Consider the quantum Hamiltonians QHi{z) defined by: det{dz — L{z)) = 
E^=o,...,nQH^iz)^l. Then: 

ykj = l,...,n andu,v eC [QHk{u),TrL^'\v)] = 0, [TrL^^\u),TrL^^\v)] = 0. (62) 



Sketch of the Proof. The proof follows immediately from the Newton identities above, theorem 



of [Ta04] (see formula[2]) on the commutativity of coefficients of det{dz — L{z)) = QHi{z)dl, and 



the binomial formula above. □ 

Corollary 1 If L{z) = Lgi^[t\{z) is given by formula fWh). thenTrl}^^ ^^1(2) are generating functions 
(in the variable z) for the elements in a commutative Bethe subalgebra in U{gln[t]). Such elements 
for k = 1, ...,n algebraically generate the Bethe subalgebra. 

Corollary 2 If Lfuii{z) is given by formula (E^j, then : TrL^^lii{z) : are generating functions (in the 
variable z ) for the elements in the center of Uc=crit{gln[t , t~'^], c). Such elements for k = 1, n alge- 
braically generate the center. (See Appendix (section [HTT]) for the definitions of Uc=crit{gln[t,t~^],c), 
the Bethe subalgebra and the normal ordering of operators :■■■:.) 



5.2.2 Application. The Newton identities for the Yangian 

We recall that if a matrix T{z) is a. nxn Lax matrix of Yangian type (see section [372]) .then the matrix 
T{z)e~^'' is a Manin matrix. We have defined (see formula HQ]) quantum powers for the Yangian type 
Lax operators as follows: 

T^p\z) =^ T{z+p)T{z + p-2) X ... X T{z + 1) = eP^^{T{z)e-^'y. (63) 

Writing the Newton identities for T^p^{z) (see formula H9l) we obtain generators in the commutative 
subalgebra of quantum integrals of motion that, to the best of our knowledge, were not known in 
the literature. 

With the notations of section IB. 1 . II page [2U] we denote by QHi{z) coefficients of the expansion of 
the "ordered" characteristic polynomial in powers of e^': 

rfef°'™(t-T(z)e-^^) = J2 t^QHk{z)e^^-''^^% (64) 

fc=0,...,n 

Note that the operator dz does not enter into the expressions QHi[z). Explicit formulas for QHi[z) 
are 

QHo{z) = i-ir qdet{T{z)), QH^-.iz) = (-1)("-) ^ grfet(T(^),,,...„,J. 

ii<...<ji 

Recall that qdet{M{z)) is defined by the formula qdet{M{z)) = Xlo-es Ilfc ^eT(k),k{z — k + l), and that 
the QHi{z)^s are the sums of principal q- minors of size i (in complete analogy with the commutative 
case, modulo substitution of minors by q- minors). 
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Theorem 6 



[QR,{z),TrT^'^{u)] = 0, [TrT^''\z),TrT^\u)] = 0, 



(65) 



kQHk{z-k + n)= J2 QHk+i{z - k + n)TrT^^{z) 

i:max(0,—k)<i<n—k 



(66) 



which is the Newton identity for the Yangian. 

Proof. The commutativity of the QHi[z) is well-known |KS81j page 114 formulas 5.6-5.7, N095 
Hence (1651) follows from (jSSD, which, in its turn, follows from theorem |H Indeed. 1351 reads: 



khk = hk+iTr{M) 

i:max(0,—k)<i<n—k 

kQHk{z)e^^-''^^^ = QHk+^{z)e^^+'-''^^^e-'^^TrT^'^{z 

i:max{0,—k)<i<n—k 



z 



ke^>^-^^^~-QHk{z-k + 



n) = e 



i:max{0,—k)<i<n—k 

('-"^^^ Y QHk+^{z-k + n)TrT^'^{z) 



i:max{0,—k)<i<n~k 

So we conclude that: kQHk{z — k + n) = Xli:max(o,-fc)<j<n-fe 

QHk+iiz-k + n)TrT^^{z). □ 



(67) 
(68) 
(69) 
(70) 



Remark 13 Let us note that for k = this follows from the Cayley-Hamilton theorem for the 
Yangian P page fIT\ by taking the trace. 



5.3 The (Quantum) MacMahon-Wronski formula 

The following identity for a s x s-matrix M over a commutative ring is called the MacMahon (and 
sometimes Wronski ( [U03][[GIQPS98] page 9)) formula: 



l/det{l -M)= Y TrS^'M, 



(71) 



n=0,...,oc 



where S'"M is n-th symmetric power of M. It can be easily verified by diagonalizing the matrix M. 
Theorem [GLZ03]^ the identity above holds true for Manin matrices (and more generally for their 
q-analogs) with the following definition of TrS'^M : 



TrS''M = l/n\ Y P™'' 

li,...ln:l<li<n 



(72) 



We remark that in this formula repeated indexes are allowed. We recall that the permanent of a 
Manin matrix was defined by formula [8] as follows: 



perraM = ■perm!''''^ M = Y^ H ^«-<^«- 

a&Sn i=l,...,n 



(73) 



This definition of traces of symmetric powers is the same as in the commutative case, with the proviso 
in mind to use row permanents. 
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see the papers Gurevich, Isaev, Ogievetsky, Pyatov, Saponov papers |GIOPS95] - [GlOPSOSj . [UOS] 
EP06] I [HL06] I [KPak06] | [FH06] | for related results 
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5.3.1 Application. Tr S^L{z), ZUcrit{gln\t,t ^],c), new Gaudin's integrals 



In this section we apply the quantum MacMahon-Wronski relation to construct further explicit gen- 
erators in the center of Uc=crit{,gln\t,t~^],c) and in the commutative Bethe subalgebra in U{gln\t\)), 
i.e., we give explicit formulas for another set of quantum conservation laws for the quantum Gaudin 



system. They are quantum analogs of TrS'^L 



classical i 



duced implicitly in [Na91 
construction. 



(the last formula - page 131, see also [U03 



The elements in the center of U {gin) intro- 
are particular cases of this 



Theorem 7 Let L{z) he a Lax matrix of Gaudin type - see definition O Let us define the el- 
ements Sn{z) as follows: write TrS'^{dz — L{z)) = ^^^g n^k,n{z)d^ and set, by definition, 
Sn{z) = (— l)"co,„(2;0 . These operators commute among themselves and with the coefficients of 
the quantum characteristic polynomial QHi{z) , defined via (det{dz — L{z)) = 'YliiQHi{z)d\). Explic- 
itly 

"ik.p and z,ueC:[Sk{z), Spin)] = [Sk{z),QHp{u)] = (74) 

Hence they provide new generators among quantum commuting integrals of motion of the Gaudin 
system. 

Corollary 3 Gonsider L{z) = Lgi^^t]{z) given by formulalWb. Then Sn{z) are generating functions 
in the variable z for the elements in the commutative Bethe subalgebra in U{gln\t\). Such elements 
for k = 1, ...,n algebraically generate the Bethe subalgebra. 

Corollary 4 Gonsider Lfuii{z) given by formula then : Sk{z) : are generating functions in the 
variable z for the elements in the center of Uc=crit{,gln\t,'t~^],c). Such elements for k = l,...,n 
algebraically generate the center. (See Appendix (section [8?T|) for definitions of Uc=crit{gln[t,t~^],c), 
the Bethe subalgebra and the normal ordering : ... :.) 

Sketch of the Proof. The proof follows from the MacMahon-Wronski relations and a theorem 
of [Ta04] (see formula [2]) on commutativity of coefficients of det{dz — L{z)). Corollary [3] follows 
immediately, corollary [3] with the help of theorem [8] (see Appendix (section l8.ll) ). □ 

Remark 14 Applying the same construction to the Yangian case i.e. TrS^{T(z)e~^'') one obtains 
expressions commuting among themselves and with coefficients of det{l — T{z)e~^^). Actually these 
expressions are quite well-known. 



6 Quantum separation of variables 

Let us briefly discuss some results and conjectures about the problem of separation of variables. 
E.Sklyanin (see surveys [Sk92] , [Sk95] ) proposed an approach which potentially should give the 
most powerful way to solve integrable systems. This program is far from being completed. The 
ultimate goal in this framework is to construct coordinates ai, Pi such that a joint eigenfunction of 
all Hamiltonians will be presented as product of functions of one variable: 

^(/3i,/32,...)=n'^'^'"*"''(^^)- 

i 

We consider this construction at the quantum level, trying to extend in this realm some ideas of 
Sklyanin andYAHH93], |DD94] . Ge95^nd others. We have checked the validity of the conjectures 



^Namely, the quantities Sn{z) are the 0-th order parts of the differential operators TrS"{dz — L{z)). 
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below for 2x2 and 3x3 cases (in particular comparing with the results of [Sk92b ). Let us remind 
that, in the classical case, the construction of separated variables for the systems we are considering 
goes, somewhat algorithmicallj0, as follows: 

Step 1 One considers, for a gin model, along the Lax matrix L{z), the matrix M = A — L{z) and 
its classical adjoint M^. 

Step 2 One takes a vector ip by means suitable linear combination of columns (or rows) of M^; in 
the simplest case, one can take ijj to be one of the columns, say the last of M^. One seeks for 
pairs {Xi,Zi) that solve 

^- = 0,z = 1, . . . ,n. 

Step 3 To actually solve this problem, one proceeds as follows. As each component ipi oi ip is a. 
polynomial of degree at most n — 1, one can form, out of ip the matrix M^, collecting the 
coefficients of the expansion of ipi in powers of A, i.e.: 

[M^]j^i = resA=o^/'iA""^"\ i = 1, . . . ,n, j = 1, . . . ,n 

Step 4 The separation coordinates are given by pairs {\i,Zi) where Zi^s are roots of Det(M^) and 
Aj are the corresponding values of Aj E§, that can be obtained, e.g., via the Cramer's rule. 
By construction , the Jacobi separation relations are the equation(s) of the spectral curve, 
Det(A-L(^) = 0. 

6.1 The Yangian case 

Let T{z) be a Lax matrix of the Yangian type (see section so (1 - e-^^T{z)]3 is a Manin 
matrix and adjoint matrix can be calculated by standard formulas (see section l^?Ti) . Consider 

/ (1 - e-^-T(z))^t°"* - (1 - e-^-T(z))fr"* \ 
{1 - e-^^T{z))'"^^°'''' = I ... ' ... ... ' . (75) 

\{l-e-'^T{z)):f'-' ... {l-e-'^T{z)rS--^ J 

Let us take an arbitrary column of this matrix, say the last column, and define Mij{z) as follows: 

;i - e-^^T{z))r;r' \ / E.=o...n-i e-'^^M,,,{z) 

;i - e-^^T(z))^^-* / V Ek=o...n-i e-'^^MnM 
Let us call "M-matrix" the matrix of these coefficients, that is : 



(76) 





M„x„(z) I I . (77) 

\ Mnfl{z) ... Mn,n-l{z) ) 

That is, Mij is a matrix collecting the coefficients of expansion in left powers of e"^'' of the elements 
of the last column of (1 - e-^^T{z)Y'^^°'''^ 

Proposition 10 For n = 2,3, M{z) is a Manin matrix. (Remark the transposition in formula [7l\)). 



^^We are herewith sweeping under the rug the problem known as "normalization of the Baker Akhiezer function" . 
^^In the quadratic i?-matrix case, actually one has to take as canonical momenta, the logarithms of these A,;, 
^^in previous sections we worked with (1 — T{z)e~^'), which is actually more well suited to the standard Yangian 
conventions, but let us switch here to (1 — e^^-r(z)) 
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6.2 The cases n = 2, 3. 



In the cases of low matrix rank, computations can be done explicitly. The following arguments hold. 

1. Define B{z) = Det^°'""''^(M(2)); then 

[B{z),B{u)]=0. (78) 

2. Consider any root P of the equation B{u) =0, (it belongs to an appropriate algebraic extension 
of the original noncommutative algebra R). Then the overdetermined system of, respectively, 
2 and 3 equations for the single variable a has a unique solution: 

,oi^)\subsU ute left z^fi 
\ Mnfi{z)\substitute left 

3. Consider all the roots (3i of the equations B{u) = 0, and the corresponding variables a^. Then: 

• The variables a^, f3i satisfy the following commutation relations: 

[ai,Pj] = -ai6ij, (80) 
[a,,a,] = 0, [(3i,(3j]=0, (81) 



Mi^n-liz)\ Substitute left 



a 



a" " \ 

n-2 



,n-l{z)\substitute left z^f3 / \ ^ / 



/ \ 





(79) 



satisfy the "quantum characteristic equation": 

{substitute left z^fi- 



\/i : detil - e ^'T{z))\ substitute left z^p,-, e-s.^a, = (82) 



• If T{z) is generic, then variables ai, f5i are "quantum coordinates" i.e. all elements of the 
algebra R can be expressed via a^, Pi and the centre (Casimirs) of the algebra R. 

The proof of these statements can be done by direct calculations; we remark that our formulas 
reproduce those of the paper [Sk92b] , 

Conjecture: It is natural to conjecture that the same hold for higher values of n, that is, that the 
M matrix is a Manin matrix, and it provides quantum separated variable for the gl{n) Yangian case. 

Remark 15 In the classical case, (and also for the Gaudin model) this solution of the Separation 
of Variables problem can be explicitly found in ,[AHH93]( The Poisson version of the conjecture above 
about Manin properties of Yangian systems, was not yet, to the best of our knowledge, considered 



in the literature; however, it can possibly be deduced from the results of [DD94j and Ge95 



Remark 16 Another open problem is a corresponding conjecture for the Gaudin case. The con- 
jecture holds true for n = 2, but seems not directly extendable to the case n > 2. 

Remark 17 We end this section with a few other remarks on separated variables. Let us consider 
the noncommutative algebra R is the algebra which enters the definition of T{z) (see section [X^ . 
i.e. by definition the matrix elements of T{z) belong to R{{z)). 

One can realize any representation of this algebra as a subspace in the space of functions F{(5i). 
Consider the joint eigen-function \l/(/3i,/32, •••) of all quantum commuting Hamiltonians 

det(l-e"^^T(^))^(A,/52,...)= ^ e~''^^hk{z)^{(3^, [32. ...) (83) 

fc=0,...,n 
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Remark that hk{z) are eigenvalues of the corresponding Hamiltonians QHi{z). 
It is easy to see from [821 that 



Vi : 



J2 e-'=^''»/ifc(A)vl/(/3i,/52,...) =0. 



(84) 



fc=0,...,n 



So the equations for functions of many variables have been separated i.e. it satisfies the system of 
equations in each variable (3i separately. 

Although from the property above one might expect that 



where each ■qj^-p<^rucie(^^^ solution of the equation ^^^e~^^''hk{z)'^^~^°''^^'"^^{z) = 0, we remark 
that formula [HI] guarantees a weaker property, namely that /32, •••) is a linear combination of 

such expressions. 

7 Concluding remarks and open problems 

We postpone further discussion to the subsequent paper(s) [CFRu] and restrict ourselves with some 
brief remarks. 

We think that it is very important to continue further the study of the notion of "quantum 
spectral curve" discovered by D. Talalaev. In particular one may hope to describe quantum action- 
angle variables with some "quantum Abel-Jacobi transform" related to Talalaev's "quantum spectral 
curve". This might also be important for the general development of noncommutative algebraic 
geometry. More generally one may hope to develop the methods used in classical integrability: Lax 
pairs, dressing transformations, tau-functions, explicit soliton and algebro-geometric solutions, etc. 
for quantum systems. 

It would be interesting to define quantum immanants ("fused T-matrices") for Uc=crit{g^n[t, t~^]i c) 
and to transfer the remarkable properties from [Okoun kov01shansky96] , {U (gin) case) to this more 
general case. In particular the task is to describe Harish-Chandra map for ZUc=crit{gin[t,t~^],c). 
(What is "loop version" of shifted Schurs ?). Moving away from the critical level towards (quantum) 
W-algebras is another task. 

We have presented results for ^f/^-Gaudin and (^/^-Heisenberg's XXX systems. We plan to extend 
our study to Q-analogues (i.e. Heisenberg's XXZ system) and super analogue in the next future, 
it would be interesting to understand a relation with more recent quantum group like structures 
( |F07] ). The challenge is to consider other systems. For example dz — Lso{n) or sp{n)-Gaudin{z) is 
not Manin matrix. However, there are evidences that there should exists an appropriate "det{dz — 
Lso{n) or sp{n)-Gaudin{z))" , siucc the existcncc of such an object is somewhat the core of the Langlands 
correspondence. The practical problem is, however, to find an explicit and viable formula for the 
determinant. In this respect we notice that, despite det^°^^"^^{dz — Lso{n) Gaudin{z)) does not depend 
on the order of columns, it does not produce commutative elements already for so(4) (A.Molev, 
A.C.)). More in general, in view of the existence of a huge number of quantum integrable systems, 
it seems to be natural to address the questions herewith studied to all of them. Namely, we point 
out the following general question: what are the conditions on a matrix with a noncommutative 
entries such that there is a natural construction of the determinant and linear algebra works? How 
efficiently can they be used in the theory of quantum integrable systems? We plan to come back to 
these and related questions in the next future. 




(85) 
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8 Appendix 

8.1 The center of Uc^crit{gln[t,t~^, c), Bethe subalgebras and normal or- 
dering 



Let us recall the main result of [CT06] and background. (See also |CMj ). 



Consider: Lf^ii{z) = ^ | | (86) 



l = — 00...00 



Here, one has to consider ei^t^ as elements of gln^®t ^gl'^\t rather than elements of gln\t-i t ^]). 
In such a case, Lfuu{z) is of Gaudin type, (see definition [5] section 13. II) . i.e., it satisfies the relations 



Via Talalaev's construction one defines the commutative subalgebra in U {gln\t\®t~^ gl^tl^'^]) ■, i-e- 
expressions QHi{z) defined by def^^dz — Lfuu{z)) = ^kd^QHk{z) are generating functions in z for 
some elements in U{gln^] © gl^\t~^]) which commute by Talalaev's theorem and hence generate 
some commutative subalgebra. 

Definition 9 This commutative subalgebra will be called "Bethe subalgebra" in U{gln[t] © 

Definition 10 Let us call by "Bethe subalgebra" in U{gln[t]) the commutative subalgebra defined in 
a similar way with the help of det'^°^{dz — Lgi^^[t]{z)) , (see IWh for Lgi^^t])- This Bethe subalgebra 
is clearly an image of the Bethe subalgebra in U{gln[t] © t~^gl'^^[t~^]) under the natural projection 

U{gl4t](Bt~'gl'^P[t-^])^U{gL[t]). 



Remark 18 The name "Bethe subalgebra" was proposed in [N095] for a related subalgebra in 
(twisted)-Yangians. (For Yangians (not twisted) the Bethe subalgebra (without a name) was defined 
in |KS81] formula 5.6 page 114.) 

Normal ordering Let us recall the standard definition of the normal ordering (see e.g. V. Kac 
|Kac97] formula 2.3.5 page 19, formula 3.1.3 page 37). 

Definition 11 Let a{z),b{z) be arbitrary formal power series with values in arbitrary associative 
algebra (in our case: a{z) = {Lfuu{z))ij,b{z) = {Lfuii{z))kh for some i,j,k,l, (see\MforLf^ii{z)). 
The normally ordered product : a{z)b{z) : is defined as follows: 

: a{z)b{z) : =^ a{z)+b{z) + b{z)a{z)., (87) 



where (e.g. V. Kac ^ Kac97^ formula 2.3.3 page 19), 

a{z)+ =^ ^ a_i_i2;* = ^ a„2:~"~^ 

j>0 n.<0 



)_ =^ Yl = E (89) 



a[z^ 

i<0 n>0 



Definition 12 (e.g. V.Kac IKacQT^ formula 3.3.1 page 42) The normally ordered product of more 
than two series a^{z),a'^{z), ...,a^{z) is defined inductively 'from right to left": 

: a\z)a\z)...a^{z) ■.=^: a\z)... : a"-^(2)a"(2) :...:. (90) 



^^by g°P is denoted Lie algebra with an opposite commutator: [gi,g2]g°p ~[5ii52]g 
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The center of Uc=crit{gln\t,t ^],c) and the "critical level" Let {gln[t,t ^],c) be the central 
extension of the Lie algebra gln[t,t~^]: 

g2t'] = [gu (72]t'+' + c{nTr{g,g2) - Tr{g^)Tr{g2))k6k,-l. (91) 

Note that for s/„ the term Tr{gi)Tr{g2) disappears and this is the standard central extension up to 
normalization. 



Fact [Ha88p, |GW89l [FF92] The center of f/c=«(^/„[t, t"^], c) ^ is trivial: Cxi, unless n ^ -1. 
For K = — 1 there exists large center, c = — 1 is called "critical level". 



Theorem 8 

J.C0I 



: det'^° {Lfuii{z) — d^) : generates the center on the critical level (92) 

i.e. define the elements Hij by : det'^°\Lfua{z) — d^) := SMmj=_oo...oo; j=o...nHijZ'^di, then Hij freely 
generate the center of Uc=-i{gln[t,t~'^],c). Here we consider Lfuii{z) (see\8^ as matrix valued gen- 
erating function for generators of Uc=-l{gln[t,t~^],c), i.e. Cmnt^ £ Uc=-l{gln[i)'t~^]-iC). 



The theorem was first proved in [CT06] , but the definition of the normal ordering given there is 



different. The coincidence with the standard normal ordering is due to |CMj . where another direct 



proof of the theorem will be given. The proof in [CT06]| is quite short, but requires various results: 



Talalaev's theorem, existence of the center, AKS- arguments, remarkable ideas from [Ryl] . [Ry3] . In 
|CMj direct proof will be given, as a corollary new proof of Talalaev's theorem will be obtained. 

8.2 The geometric Langlands correspondence over C 

The geometric Langlands correspondence over C [BP] (excellent survey [Fr95]') is related to the con- 
structions above as follows. (See JCT06]^ which originates from _ [FFR94] , [GNR99] ^ and is intimately 
related to Sklyanin's separation of variables). As in number theory there are two versions: local and 
global. 

8.2.1 The local correspondence: 

" Automorphic side": Irreducible Representation V of {gln[t,t~^],c)c=crit 

I Local Langlands correspondence over C \ 
"Galois side": differential operator of order n: J2k=o n^z^k{z) (93) 

In number theory the local correspondence is a bijection between Irreps{GLn{Qp)) and n-dimensional 
representations of the Galois group Gal{Qp, Qp), where Qp is the field of p-adic numbers. The stan- 
dard analogies: Qp is analogous to C((t)), irreps of GL„(Qp) to irreps of {gln[t, c), representations 
of the Galois group of Qp to differentials operators (since Galois group is similar to the fundamen- 
tal group (i.e. it is unramified part of the Galois group, when both are defined) and a differential 
operators defines a representation of a fundamental group via monodromy of its solution). 



Construction CT06 : Consider irrep V of {gln[t,t~^],c)c=crit, by "Schur's lemma" center of 
Uc=crit{gln[t,t~^], c) acts by scalar operators on any irrep, hence one obtains a character xv of the 
center xv '■ ZUc=crit{gln[t,t~^], c) — > C. Now the desired differential operator is given by Talalaev's 
formula: Xv{- {.det^°\dz — Lfuii{z)) :). (See formula [86] for Lf„u{z) and[92l) Summary: 

V G Irreps{gln[t,t-\c) ^ Xv : ZU,=,U9ln[t,t-\c) ^ C ^ Xv{- {def°\d, - Lfuii{z)) : ). (94) 
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more precisely one needs "local completion" of Uc=K{gUi[t^ t c), see [FFR94]| 
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The non-triviality is hidden in Talalaev's formula and theorem [HI The correspondence is important 
in integrability since analysis of {det'^°\dz — L{z))Q{z) = allows to write Bethe equations for 
the spectrum and defines eigenf unctions in separated variables. It also satisfies many geometric, 
categorical and "stringy" properties, but this is out of our scope. 

8.2.2 The global correspondence 

Consider an algebraic curve S, denote by BuHqi^^ the moduli stack of GLn vector bundles on S. 
The original version of the correspondence [BP] is the following (later it was extended): 

" Automorphic side" : Hitchin's D-modules on Bun^j^^ 

I Global Langlands correspondence over C \ 
"Galois side" : GLn — opers on S 

i.e. roughly speaking differential operators of order n on S (95) 

The construction (not the proof) of the global case is actually an immediate corollary of the local 
one ( |BDj ). It goes as follows. 

Denote by HH commutative algebra of differential operators on Bun^^^ which are quantiza- 



[Hi^ 


^7], 


C. 


By 



by characters xd '■ HH C By nowadays standard arguments (i.e. using representation 
Bun^^^ = Gin\GLn{{z)) / Gout) one sees that the commutative algebra HH is actually a factor 
algebra of ZUc=crit{gln[t,t~^],c). So one gets xd '■ ZUc=crit{gln[t,t~^], c) C and by the local cor- 
respondence one gets the differential operator Dz. One needs to check that it is actually globally 
defined object on E |BD] . (It is quantization of a globally defined object - spectral curve, so hopefully 
in future it would be derived from general facts). □ 

Remark 19 For S = CP^ with the "marked points" ("tame ramification") the Hitchin system 
(H.s.) is precisely the standard Gaudin system with Lax matrix [TTl for marked points with multiplic- 



ities ("wild ramification") - [T8b (see |Mar94j . |[ER95 , Ne95 ). For S = CP^ with nodes and cusps 
see CT03A , more general singularities [CT03B] , H.s. can be seen as Knizhnnik-Zamolodchikov 



equation degenerated to the critical level. 



8.2.3 Remarks 

Let us emphasize the two features of the integrability point of view on the correspondence. 

• The correspondence has nothing specific to Hitchin's integrable system, it should work for 
an arbitrary integrable system which has a spectral curve. The analog of G-oper should be 
"quantum spectral curve" defined by Talalaev's type formula det{X — L{z)). (For example 
in Heisenberg XXX-magnet one should write Det{l — e^^-T(z)), in super version: Ber{l — 
e-^^T(.)), ...) 

• The geometric Langlands correspondence over C is just an appropriately understood quanti- 
zation of the basic correspondence in integrability: Liouville tori of classical integrable system 
correspond to some spectral curves. Quantization sends Liouville tori [H^ = Aj) to Hitchin's D- 
module (essentially one needs to construct Hi, the D-modules are generated by Hi — x{Hi), for 
an arbitrary x '■ Hi ^ C) and (amusingly) the spectral curve itself should be quantized by Ta- 
lalaev's type formula det{X — L{z)). The quantum spectral curve is a differential (difference - de- 
pending on system) operator in one variable. The corresponding equation det{\— L{z))Q{z) = 
generalizes Baxter's T-Q relation. It is a key ingredient in a solution of quantum systems. 

^^actually twistcd-difFcrcntial operators 
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The geometric Langlands correspondence is rapidly developing in various directions, much work 
should be done, one may hope for exciting achievements. 



8.3 Capelli and Baxter type conjectures 

Let us announce some conjectures on Talalaev's quantum spectral curve, the Capelli identities and 
the Baxter type conjectures. We plan to discuss this in subsequent papers. 



8.3.1 det(dz — L{z)) and det{dz + L{z)) are conjugate 

Let L{z), T{z) be Lax matrices of the Gaudin and Yangian types (see definitions [5], page [TO] and [6| 
page [12] respectively). Under the formal conjugation D D* we mean the following: (^^ fi{z)d\)* = 
"^ii—dzY fi{z) i.e. the conjugation does not act on coefficients /j. 

Proposition CI The differential (difference) operators det''°\dz — L{z)) and (—l)"'det^°'^{dz + 
L{z)) (respectively det^°\l — T{z)e~^^) and def°'^{l—e^^T{z))) are formally conjugate to each other. 

Sketch of the Proof. The Yangian case follows from the known result, Gaudin case follows 
from the previous one by degeneration. Indeed, obviously: 



det'-\l-T{z)e-'')= J2 E ?^et-'(T(%,....,J(-e)- 

k=0,...,nji<...<jk 

def°'"{l-e^'T{z))= {-ef^^qdef'""{T{z)j,^,„ 



(96) 
(97) 



fc=0,...,n ji<...<jfc 



where qdet""^ {M (z)) 

So-e5„ nfc=i,...,' 
qdet'^\M{z)) = qdet 



n 



fc=i, 



k 



qdet''°'"{M{z)) 



n -\ 
{M{z)) 



k - 1) 

[see e.g. formulas 2.24, 2.25 page 10[Mo02 ). Also recalling that 



For the Yangian type matrix M{z) it is true that: 



any submatrix of the Yangian type matrix is again Yangian type matrix we come to the desired 
conclusion for the Yangian case. 

Now we follow the ideas of D. Talalaev | [Ta04] ': one can represent T{z) = l + hL(z) + h'^ (...), hence: 



1 - e^^^T{z) = -h{dz + L{z)) + h^{...). So det'°^{l - T{z)e-^~~) = {-hYdet''°\dz + L{z)) + h^'+^i...). 
Similar with 1 - T{z)e-^^'-. We proved above that det^°\l - T{z)e-^^') an d &r°" '(l - e^^-T{z)) 
are formally conjugate (it was proved for /i = 1, but one can rescale with /i[Ta04 ). So their first 



coefficient in expansion in h are formally conjugate, that is the desired result for the Gaudin case. □ 



8.3.2 Capelli type explicit formulas for det(dz ± L(z)) 

Conjecture C2 The Capelli type explicit formulas. Assume thatnxn Gaudin type Lax matrix L{z) 
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has only simple poles at some points z = Zi,i = 1, k (e.g. formulas [7^ [77| j, then: 

det^°\d, - L{z)) = {def'"^{d, + L{z))y~^^y conjugate ^ 

(98) 

= E (-1^ E ciet™'(M+ J] -^dzag{z-l,z-2,...,0})d:-' = 

principal i x z i;=i,.--,fc 
submatrices M of L{z) 

(99) 

= E E cier'(M+ 5^ ^_rf,a^7{0,-l,...,-z + l}) = 

principal i x -i j;=i,.-,fc 
submatrices M of L[z) 

(100) 

A^'ow row-dets: 

= E E c^er°'"(M+ J] ^rf^a^7{0,l,...,z-l})ar = 

principal i x i v=i,...,k 
submatrices M of L{z) 

(101) 

= E (ier°"'(M+ ^ _l_rf^a^?{-z + l,...,-l,0}). 

principal i x z i)=i,...,fc 
submatrices M of L{z) 

(102) 

Analogously let us formulate conjecture on the explicit expressions for the Nazarov type elements. 
Here determinants are substituted by permanents. 

Tr(^-(9.+L(.)))=i y: eh (^--.)+ E 

l<ii<"-,l<i2<",---,l<ifc<"-o-6Sp r=l,...,p v=l,...,k 



8.3.3 General Capelli identities 

Let us denote by L{z) Lax matrix of the Gaudin type given by the formula [39l T^-^^(z) is a Lax 
matrix given by formula [241 recall that the linear map "Wick" was defined by formula HU Consider 
Hkj defined as det{X — L{z)) = J2k=o n-j=-n oo -^'^^"'-^fcJ- function made of Hkj is called 

classical integral of motion for the Gaudin system. Similar for the T^^^{z). 
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Conjecture C3 For all classical integrals of motion Hi,H2 of the Gaudin and XXX systems: 

Wick{Tr{X - L'^^^^'^'^^z))'') = Tr{d, - (105) 

Wick{TrS^{\ - L''^'^'^""'\z))) = TrS\d, - (106) 

the Yangian type systems: 

[WickiH^^''''""'^ ^^^{p, q)), WickiH^^"''''"^ ^^^{p. q))] = 0, (107) 

Wick{det{l - e-^T"''^^^^™^ ^^^{z))) = det{l - e-^^T^"""*"™ ^^^{z)), (108) 

Wick{Tr{l - e-^r^"'''""^ ^^^{z)f) = Tr(l - g-^'-T"'"^"*"™ ^^^(^))^ (109) 

Wick{TrS\l - e-^r^''''''"^ ^^^(^))) = TrS\l - e-9-r^"™*«™ ^^^(^)), (HO) 

V Young tableux : Wtcki{PrYoung taUeu.T '''''' ''''''\z)T '''''' '''''\z)....)) = 

{Prvoung tableu.T'^'''''''''^ '''''' {z + Ci ) T ^^^(z + C^)...). (Ill) 

- No "Tr" generalizes [ Ok96B]\ if one puts "Tr" - weaker [ Ok96]\ 



Remark 20 The expressions (PrYoung diagramT ^^ ^{z + ci)T9"""*""" ^^^{z + ca)....) are 
called "fused T-matrices" in physics literature. See Ok96 , [[Na96] for the concise exposition and 



deep results for the particular case of U{gln), the elements in ZU{gln) obtained in this way are called 
"quantum immanants" in loc.cit. It seems there is no concise mathematical exposition (and probably 
understanding) of general "fusion". PrYoung taUeux G Main ® ... ® Main is the Young projector on 
a irreducible gin representation corresponding to a Young tableux, numbers Cj = rowi — columni 
for the cell with a content i. The conjecture above generalizes the results from |[Ok96]^ .[Ok96B]j 
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8.3.4 Generalized Baxter type conjectures 

Conjecture C4 Let (vr, y) be a finite- dimensional representation of gl[t] or Yangian, then all solu- 
tions of the equations below are single-valued functions (products of rational function over exponen- 
tials): 

n{Tr{d, - L''"'^"*""^(2))'=)Qi(2) = 0, (112) 

n{TrS^{d, - L«"""*"™(z)))Q2(;2) = 0, (113) 

7r(Tr(l - e-^'T^''''3'°'''{z)f)Q^{z) = 0, (114) 

7r(rr^'=(l - e-^^T^'^"^''^"(2)))Q4(^) = 0, (115) 
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and more generally V Young tableux: 

1 2 

n{Wzck{Tr{PrYoung taUeu. (l " A-1T^'™'(^)) (l - a-1t^'™'(z)) ....)))Q,{z) = 0, (116) 

1 2 

7r(W^^cA;(Tr(Pry„,„, (a - L'='™'(z)) (a - L^'™'(^)) ....)))Q^{z) = 0. (117) 

This conjecture will hopefully lead to new forms of Bethe equations for the spectrum, similar to 
the considerations with the standard Baxter equations. 
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